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, The M— fold degeneracy of the ground state of a quasi-one dimensional system allows it to support 

' topological excitations such as solitons. We study the combined effects of Coulomb interactions and 

\ confinement due to interchain coupling on the statistics and thermodynamics of such defects. We 

concentrate on a two-dimensional case with M — 2 which may correspond to monolayers of doped 
polyacetylene-type polymers or other charge density waves and to junctions in field effect experiments 
for equivalent materials. Combining analytical and numerical methods, the theory is developed by 
a mapping to the 2— dimensional Ising model with long range 4-spin interactions. We obtain the 
phase diagram depending on the ratio of Coulomb interaction and interchain coupling with respect 
to temperature. The latter exhibits deconfined phases for liquids and Wigner crystals of kinks, and 
confined ones for bikinks. Also we find aggregated phases with either infinite domain walls of kinks 
or finite rods of bikinks. Roughening effects due to both temperature and Coulomb repulsion are 



I. INTRODUCTION 



■ observed. Applications may concern the melting of stripes in doped correlated materials. 

B 

> , 

Charged interfaces, from their microscopic origin to their macroscopic statistical properties, appear naturally in the 
study of doped correlated insulators. The latter are characterized by discrete symmetry broken phases with M-fold 
degenerate ground states, M > 1. The interfaces, which are topological configurations whether charged or not, 
provide a connection between the different ground states. In this framework, quasi-one-dimensional systems, e.g. 
systems of weakfef bound chains, exhibit them at a microscopic level via solitons which emerge as single particles, 
the tp— particlesau. This name is related to the fact that each particle provides a phase jump tp — 2n/M, connecting 
two different but equivalent ground states, within a chain. Their presence therefore breaks the coherence between 
neighboring chains which must have the same phase modulo 2tt. Actually they may exist in a broader range of 
materials even though their individual detection might be difficult. In this respect, -a. defects, i.e. *-M — 2, may be 
interpreted as the building blocks of the experimentally observed-stripes in manganitea and cuprateutj oxides as well 
as of the domain walls observed in uniaxial ferroelectric systernsB. The latter references,U andta, should also provide 
informations about the charge carried by kinks via measurements of the tilt angle of the stripes and the dipole 
| momentum, respectively. 

£N) ■ In one dimension (D) or above the transition temperature, uncharged ^—particles arc allowed to exist as a gas of 
elementary excitations. They can be defined as strings of reversed spins which, in D > 1, cost an energy propor- 
tional to their length because of interchain coupling. Below the transition temperature T c , in the ordered broken 
symmetry phase, a single soliton would cost infinite energy. The condition for their preservation lies in the concept of 
confinement, which limits the string's length leading to finite energy excitations. T c is thus called the confinement 
- deconfinement transition temperature. Below T c , a constant force is originated by confinement which binds the 
kinks into trains of M particles, to recover the total increment of 2tt, within each chain. Going to zero temperature, 
, confinement leads to the aggregation of the trains of solitons belonging to different chains. Thus, the microscopic 
defects become mesoscopic complexes which span the whole sample at T — and can be viewed as domain lines or 
Q , walls whether the system is, respectively, 2D or 3DB. 

These qualitative arguments strictly apply when solitons are uncharged and must be extended to the case where 
, long ranged Coulomb interactions are present. Indeed, doping correlated insulators leads to the formation of charged 
topological defectsO. Also, field effect experiments^ provide junctions of mobile charged kinks interacting with their 
rS images. 

The present work deals exclusively with a 2-fold degenerate ground state which corresponds, for quasi-one-dimensional 



systems, to the well known case of polyacetylene, cf.oEInij. In such (M = 2) systems the elementary defects are the 
±7r-solitons. The kink (+) shifts the phase of the CDW from to tt while the anti-kink (— ) shifts it from n to 27T. 
Both amplitudc-solitons have the same electric charge. The anti-kink of opposite charge, corresponding to the tt to 
shift, is not thermally activated. For the sake of clarity we thus attribute to the — 7r-soliton the same graphical 
representation as the former in Fig. [j]. The trains of kinks are then simply pairs of kink— anti-kinks or "bikinks" and 
the domain lines perpendicular to the chains are formed by such bikinks. We shall examine the statistical properties 
of kinks under the combined action of the attractive confinement force and the repulsive 3D Coulomb force with the 
help of both analytical and numerical methods. 
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The paper is organized as follows. In chapter || we present the general model of confinement based on the mapping 
of a system of interacting kinks to the anisotropic Ising model. In chapter II] the phase diagram of the 2D system, 
with competing Coulomb interactions and interchain coupling, is given. Analytically, the results are obtained by 
various estimates in the case of strong correlations and with the help of a phenomenological model, equivalent to the 
Ising model below T c , which has been extended to include weak Coulomb interactions. The thermal roughening of 
uncharged domain lines, as well as the T — Coulomb roughening, arc also considered. In chapter [V, numerical 
simulations are presented. Various regimes of the estimated phase diagram are discussed. The impact of discreteness 
imposed by numerical simulations is also considered. Details of calculations are given in the Appendix. 



II. THE GENERAL MODEL OF CONFINEMENT 



We shall follow the method oiEl, where the statistics of uncharged kinks in a system of weakly bound chains was 
studied by mapping to the Ising model. 

The original Hamiltonian describing the kinks reads 

H = H - [xN s = j dx{~ Vr] a (x)T]ff(x) + ^(£ s - fJ, s )c a {x)}, (1) 

where N s is the total number of kinks, E s is the energy of an isolated kink, fi s and c a (x) are, respectively, their 
chemical potential and local concentration. V is the interchain coupling constant, chosen to be positive. n a (x) is the 
order parameter describing a kink along chain a. 

On the atomic scale ao x , the kink extends over a core length £o > ao x , e.g. £o ~ 7acte f° r polyacetylene, around 
position x. For this system, the microscopic, domain-wall like, excitation takes the values n = ±1 between the kinks 
and far from their cores. 

The mapping to the two-dimensional Ising model, whose Z2 symmetry is clearly compatible with the 2— fold degeneracy 
of our system, is then based on a coarse graining procedure by introducing a new length scale a x > £0 > a 0x , which 
corresponds to the optimal width of a pair of kinks within a chain under the combined action of interchain interaction 
and exponential intrachain repulsion between them, cf. Fig. |l|. We then discretize the system by and 
y — ► y a = aa y where a y is the interchain distance. On this scale the order parameter can be restricted to the values 
±1 so it is written as an Ising spin, S n ^ a . The kink density, corresponding to the number of kinks per site, can also be 
expressed in terms of the Ising spins in the following way. The kink, on chain a for example, is defined between two 
consecutive opposite spins, say on sites n and n + 1, cf. Fig. [j]. The dimensionless density of solitons p nyCC is related 
to their dimensional concentration c a (x), cf. (|l|), by 



In terms of spins these considerations yield 



C-a{x^CL x — Pn,a- (2) 



<SVi,a<SVi+l,a)- (3) 



We see that p n ,a is unity when a kink is present and zero otherwise. 

Introducing these new variables in the first Hamiltonian (|l]) and including Coulomb interactions we recover the 
two-dimensional Ising model 

H = Hq + H c = —J± ^ S n , a Sn,l3 — J\\ ^(Sn,aS n +l,a — 1) + Hcoul, (4) 
<a,0>n a.n 



where 



J ± =Va x =K ± T J\\ = l(E a -n a ) = K\\T. (5) 
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The first term in (||) describes the interchain interaction, responsible for the confinement of the kinks, whose strength 
is given by the dimensionless coupling constant K±. In the second term, K\\ is the effective chemical potential of a 
kink, which has been included as a coupling constant. We emphasize the fact that the energy of a free soliton E s is 
always larger than the temperature, so that none of them is thermally activated. We only consider redistributions of 
the kinks, no matter how they have been created (doping, charge transfer, incommensurabilities), under the combined 
action of confinement and Coulomb interactions. The latter has been introduced in (^]) through 

H C oui = Ye 2^ ~\T -r> 7\ ' () 

n,m;a.p 



where 'r'n.a = {na x , aa v ), v = c a (x)a x is the average dimensionless density of kinks, see (^) below, and e the dielectric 
constant of the isotropic, neutral media in which the plane is embedded. The Coulomb interaction as written above 
takes into account the presence of a homogeneous, negatively charged, background of density v which ensures the 
electroneutrality of the whole system. It gives the total electrostatic energy of the kinks and the background they are 
interacting with. The case of dipoles, produced in field effect junctions will be considered later. 

As we are working in the grand-canonical ensemble, the key for finding the thermodynamic properties of the system, 
is to determine the chemical potential of the kinks via the average density per site of kinks: 

»=5<'-<s~w »-i( i+ ?J^)— < 7 > 

where the averaging < . >_a | nd the free energy // are defined for the ensemble of Ising spins and Q s is the grand 
potential of the solitons, cf.EJ. 
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FIG. 1. Figure a) displays a tv— soliton (solid line) along a chain (broken line). The arrow indicates the width £o of the defect. 
Figure b) displays the soliton and the antisoliton which correspond to strings of inverted Ising spins after coarse 
graining. The anti-soliton corresponding to a shift of the phase from tt to 2tt has the same charge as the soliton. 
The other anti-soliton (dotted line) has an opposite charge and is not thermally activated. 
Figure c) displays a soliton and an antisoliton of the same charge. For simplicity the ir to 2ir kink has been 
represented as the tt to one which does not exist in our case. Neighboring chains (broken lines) are supposed to be 
free from defects. Due to interchain interactions (bold lines) the energy of the bikink increases linearly with its size. 
As a consequence, a constant confinement force binds the solitons. Exponential core repulsion prevents them from 

annihilating. 
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III. THE TWO-DIMENSIONAL SYSTEM 



A. The dimensionless parameters of the phase diagram 



Without Coulomb interaction, the Ising model mapping allows for a complete analysis of the 2D case. The aggregation 
of kinks is controlled by the temperature or the reduced temperature used in the phase diagram, 

If their density is weak enough, which we will suppose in the following, below the confinement transition temperature 

T c «^ , c «i (8) 
v v 

the growing confinement energy binds the kinks in pairs with remaining energy Jj_ . When the bikinks on neighboring 
chains stick together below the crossover temperature, 

T ° ~ i — 7^~\ > ®° ~ i — TTTTj (9) 
log(2z/) log(2i/) 

the energy J± is released which results in an adhesion process leading to growing rods. The latter are double lines 
formed by aggregated bikinks, perpendicular to the chains. Their growth when temperature is further decreased can 
be seen from the density of solitons and the average rod length l± 

exp(-4Xj_) 



, lx = ^- - OyV2u exp(2^ x ). (10) 

These expressions, obtained from the exact solution of the non-interacting Ising modeJll show that as T — > the 
length of the aggregates grows exponentially with temperature until they cross the whole sample more like a system 
of free domain lines. 

On the other hand, the strength of the Coulomb interaction can be characterized by 

r = — Wy = ^^ (ii) 

J± ea y 

where s, = a x a y is the unit area of the coarse grained lattice and z < 1 takes into account possible fractional 
charges^ which may be measured in ferroelectric experimental! Parameters Jiand w y give, respectively, the scale 
of the confinement energy and of the Coulomb energy perpendicular to the chains. They are invariant under coarse 
graining, thus T is also invariant. Unless e is very large the w y scale is larger than other ones. For example, it is much 
larger than the energy per particle of a Wigner crystal, which is of the order of e 2 /eR where R ^> a x . y is the mean 
distance between particles of the Wigner crystal. It is also important with respect to the confinement. In Rydbergs, 
lRy = e 2 /eao, where ao is Bohr's radius, we have w y = za^Ry / ea y . In the same units, Jj_ ~ vT c w vl0~ 3 Ry if we 
choose T c as 100K. Thus for a y = 5ao and v — 1% we have 

r = 3i w Z 10 2 /e. (12) 
J± 

Taking z — 1, the dielectric constant needs to be of the order of 10 2 for the confinement energy to be of the order of the 
Coulomb energy. This is large but may not be excluded for e ~ (u> p / A) 2 , where the plasma frequency to p « 10 _1 — 10 
eV and the gap A w 10~ 2 - 10° eV. 

At low temperatures T controls the aggregation of the defects in a way similar to 9 at high temperatures. This follows 
from the derivation of the optimal width x op t of a bikink under the combined action of Coulomb force e 2 /ex 2 and the 
confinement force —V: 

x op t = vsT (13) 

This relation holds in all bikink phases whether th ey a re liquid or of the Wigner crystal type. In both phases the 
average distance between bikinks is given by R — \fsjv. As long as x op t -C R the elongation x op t can be neglected so 
that, in zero approximation, the bikink can be considered as a well defined particle of charge 2e. On the other hand 
when x op t is of the order of R, bikinks deconfine and the system consists of interacting solitons. These remarks are 
weakly modified at non zero temperature as long as 9 < 1. This follows from the fact that thermal fluctuations shift 
the mean width of the bikink to x ~ x op t + a K r 1 / 4 9 1 / 2 . Thermal fluctuations can be neglected as soon as T ^> O 2 . 



Unless e is large, cf. (12), this condition is well satisfied below 6q. 
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B. Wigner crystals of kinks and bikinks 



Below 80, when the condition T ^> 8 2 is satisfied, bikinks order in the form of a Wigner crystal. The ratio 
Xtot/R — y/Tv is another control parameter which defines a crossover to a deconfined phase consisting of a Wigner 
crystal of solitons. The Lindeman criterion specifies that the crossover occurs when the mean square deviation of the 
width x op t becomes of the order of a fraction of R. With the small Lindeman number 7 the transition is thus expected 
at 

r c ~ 7/", (14) 

which scales as 9 C . 

This deconfinement transition with respect to the bikinks can be easily derived from the Wigner crystal phase of 
solitons. The energy per particle of the Wigner crystal of solitons is —e 2 /eR and the confinement energy ~ J±R/a y . 
A confinement crossover, with respect to the solitons, takes place when both energies are of the same order. Taking 
into account the Lindeman number, we recover (|l4|). 
These results are summarized in the phase diagram of figure H. 



C. The interacting rod phase 

The Coulomb interactions are expected to be especially important for aggregated phases, e.g. rods and lines, 
because of the size-dependent charging energy cost of the aggregation, 

W \~o 1 ,min(na y ,l D ) 

S(n) w 2w y n\og{ y - ), (15) 

a x 

where n is the number of solitons belonging to each of the two lines formed by the rods and l]j is the Debye screening 
length. The latter is determined by the density response function dis/dfi s which, for the 2D plane embedded in 3D 
space, is 

1 2ne 2 dv , x 

Id es dfi s 

We shall consider the case of relatively weak Coulomb interactions which does not prevent the actual existence of 
the aggregates while affecting their statistical properties. The required condition is that r<l. Our aim is then to 
formulate the general expressions which yield the density of kinks and their mean length. This can be done with the 
help of a phcnomcnological model which has been shown to adequately reproduce the results of the exact solutionEj. 
The rods can fuse into longer ones or disintegrate into shorter ones, equilibrating their chemical potentials. The main 
effect of the Coulomb interaction is to shift the energy E(n) of each aggregate of length na y by its self-energy jl5| ) 
while preserving their shapes. The free energy of the system is then 

00 

/rod - ( c ( n ) E (n) + Tc ( n ) log(c(n)/e)) , (17) 

where c(n) is the density of rods of length na y , the second term is the entropy of the Boltzmann-like rod gas, and 

E(n) = 4J X + 2E s n + T,(n) (18) 

is the energy of a rod of length na y . The first term in ( |l8| ) corresponds to the confinement energy paid by the tips 
of the rod. In the second term we have taken into account the fact that there are two kinks per unit length of the 
aggregate. 

Minimizing the free-energy with the constraint that the total number of particles is constant, 







frod f-^s ^ 

2nc(n) , (19) 



dc{n) 

the density of kinks is then given by 
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v = ^ 2nc(n) = e - iK± ^ 2n exp {-4nK\\ - /3£(ra)) , (20) 

n— 1 n— 1 

where is related to the chemical potential of the solitons fj, s by (|J). Thus, expression ( |20| ) allows one to fix fj, s . 
Finally, the average length of the domain line reads 

oo 

l± = a y — c=^^c(n). (21) 

n=l 

From (|2C| ) and (E^) with £(n) = we recover the results in (|l(]). As will be shown in subsequent paragraphs, 
screening by normal carriers, cf. Id, is important in the present regime. Self-screening in the ensemble of solitons 
is also considered. It plays a dominant role at T = 0, especially for stronger Coulomb interactions corresponding to 

r> 1. 



1 . The case of efficient screening 

The efficiency of the bulk screening is manifested by the fact that the Debye length lr> is shorter than the average 
l eng th of the domain line l±. Considering a distribution of rods with lengths na y > Id, we arrive, with the help of 
( |l5| ) and (|2(]), at the expressions given in ( |To| ) with K\\ — > K«, where 



K \\= K \\ + ^ w v l0 S( l f) C- ) -' 



is the effective chemical potential shifted by the Coulomb self-energy. We can now check our original hypothesis 
i_L > Id- We have with the help of ( |l6| ) 

x _ 2ire 2 exp(-4^ ± ) 1 
D ~ es T (2^,p3- [ 1 

l±_ _ 2ne 2 exp(4K ± ) 2 

which grows as T —> 0, i.e. 1±/Id —> oo. Therefore, thanks to the presence of an effective bulk screening, the 
interacting system can be mapped to the non-interacting one, in the limit Kj_ 3> 1, proving the growing macroscopic 
aggregation of the defects. 



Thus 



2. The case of weak screening 

More surprising is the case of weak screening which corresponds to l± < Id- Considering a distribution of rods with 
lengths na y < Id, ( p0[ ) then becomes 

u^e- iK ^ I" dn2nexp(-4nK n -2n^log( 1 ^)) (24) 
Jo \ " T a x J 

where the second term in the exponent comes from the self-energy of a rod of length n. Due to the non-linearity of the 
self-energy the chemical potential cannot simply be shifted as in the previous case. Moreover, the Coulomb energy 
makes the sum converge even for a zero value of the chemical potential. In order to balance the Coulomb interaction 
and thus keep a given non-zero y, we must allow K» to acquire negative values in (p4|). To show this, we notice that 
the integral in ( p4[ ) can be evaluated with the help of a saddle point approximation, which is valid for large, negative 
if I] . At the saddle point the argument of the exponential in (|24|) has a sharp maximum for 

^ exp (4r) (25) 



G 



and ( p^ ) becomes 




uexp{4K x )^2J— ^exp(^Z sp ) . (26) 



From ( |16| ) and (24) at the saddle point, the Debye length reads 

1 w. 



, - T vls P - (27) 
Id 1 a x 



As we work at constant v, l sp is given by ( |2q ) 

Jj_ T ( v w. 



I ^4^ + — log 7 ^= ■ (28) 

Due to the weak Coulomb interaction T 1 and the sharpness of the saddle point ensured by w y ^> T, the second 
term in ( |2§| ) is a correction with respect to the first. The second term is responsible for the increase of l sp with 
decreasing temperature as long as 

2 2 

T>—f=T l . (29) 



Hence, in this regime, l sp « AJ±/w y and (pTj) yields 



l sp 16J : 



2 



-za (30) 



Id Tw y 

When l sp reaches Id a crossover from the actual process of growth, involving negative Ku , to the previous one will 
take place at T cr defined, with the help of (^), by 

^ « IS" 2 -, (31) 

i c Wy 

where (ph has been used. From (E9) and (p|) we have 



T rr 16 2 / J x 3 



— »1- (32) 



where we have also used the fact that we are in the weak Coulomb regime. Thus (|32j) indicates that the well screened 
regime ( l sp > Id) is reached before temperature is low enough to disfavor the growth of the rods, i.e. (^9| ) is satisfied. 
We summarize the case of weak screening as follows. 

Without Coulomb interactions, the rods started to form at T < T ps J±/\ log f| and then grew exponentially with 
an activation energy 4J_l; as temperature approached zero, the chemical potential approached zero exponentially, cf. 

In the case of weak Coulomb interaction (w y < J±_), at low temperatures (T <C w y ) the behavior of the system is 
very different. Remarkably, negative Ku corresponds to the Ising model with antifcrromagnetic coupling along the 
chains. Nevertheless, 4— spin Coulomb interactions (^) keep the system in the ferromagnetic state so as to preserve 
electroneutrality. The average length of the aggregates l sp is approximately Jj_/w y with weak deviations linear in 
temperature up to some logarithmic corrections, cf. (|2^). These deviations lead to an increase of l sp when T decreases. 
The chemical potential T\Kn \ moves down within the negative regime but logarithmically slow as can be seen from (|2^). 
Meanwhile the compressibility « dv/dfj, s increases so that the screening length decreases linearly with temperature, 



cf. (27). At T cr the Coulomb interaction becomes screened. The system then returns to the exponential growth of lj 



as for the noninteracting system. 



3. Application to field effect junctions 

To finish with this section we mention that the logarithmically growing self-energy is peculiar to a doped system. 
In the case of a layer produced by a field effect experiment, dipoles order in lines of length a y n and the associated 
energy is 
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/ \ (ed) 2 ( 1 1 \ 

where ci is the distance between a charge and its image. As we are interested in large n, this energy becomes linear 
in n. Thus this case is similar to the well screened case. 



D. Thermal and Coulomb roughening of domain lines 



1. Thermal roughening of lines 



The effect of thermal fluctuations on the aggregates of the non-interacting_system can be derived from the solid- 
on-solid (SOS) model which is well known in the field of random interfaces, cf.ta for more details. In this model a line 
made of one topological charge, i.e. a domain line, is described by integer-valued height variables Xk (k = 1, ...,l±) 
which define the position of kinks along the chains with respect to— their T = position. The energy of an SOS 
interface at T — > 0, i.e. K± — > oo and l± — > oo, can then be taken asE_a 



pH S os = 2K ± J2\xk ~ x k -i\- (33) 
The difference correlation function is thenO 



fc=i 



Expression (53) shows that the correlation function diverges in the limit \k — j\ — > oo. This is a signature of the 
roughening of the domain line at any non zero temperature. This roughening leads to collisions between different lines 
separated by a distance a x /v. Taking into account the fluctuations of one line, neighboring lines remaining straight for 
simplicity, a collision takes place once the mean fluctuation in x— direction is of the order of a x jv. Collisions of lines 
can be equivalently interpreted as disintegration of rods as can be seen from (j34|). With < {xk — Xj) 2 >sos~ ^/v 2 , 
the temperature at which rods of length lj_ w \k — j\ are formed is given by 

^ (35) 



k>g('±- 2 ) 



Equation ( pq ) shows that rods reduce to bikinks (l± ~ 1) when T > To in agreement with what has been said in the 
preceding paragraphs. 



2. Coulomb roughening of lines 



In previous sections we have shown that the growing of the aggregates as we approach zero temperature is mainly 
due to an effective bulk screening. However, at T — all kinks have aggregated into domain lines which cross the 
whole sample. Thus, no screening charge remains and no confinement energy, due to interchain interactions, is paid 
by the lines. The latter implies that two domain lines, each made of one type of topological charge, will be far apart. 
On the other hand, because there is no more screening charge, the competition between Coulomb interaction and 
confinement leads to roughened lines, as shown below. It is a signature of the self-screening of the domain lines. This 
roughening may be considered separately from the thermal roughening we discussed previously. 

First, we give some simple estimates concerning the form of the line. The latter corresponds to a linear distribution 
of charges and it's equilibrium shape is determined by allowing one of it's constituent kinks to find it's equilibrium 
position under the combined action of the Coulomb repulsion and the attractive confinement force. The Coulomb 
force F cou i = w y /x originates from the infinite line and is balanced by the confinement force F con f = —V. The 
optimal equilibrium position of a kink emitted from the line is 

x ~ ^r-a x = Ta x . (36) 
J± 

For r> 1, the linear distribution acquires a width io > a^. This finite width is a manifestation of the self-screening 
of the rods which leads to their roughening. 
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A particular distribution of kinks within the roughened domain line requires a more elaborate analysis. We consider 
a general distribution characterized by the variational parameters xq and yo- The latter correspond to the periodicity 
of the distribution while the former is the modulation amplitude. They are found by minimizing the total energy per 
chain, E to t, of the distribution. The latter is 

Etot = E cou i + E con f (37) 

with the Coulomb energy per chain 



w y f dydy' 
L 



a l y \y - y>y + ( X (y) - X {y')f 



(38) 



where L is the size of the system in the y direction. The confinement energy per chain is 

a y J± f Vo dx 
E conf = — / — dy. 
yo a,x Jo dy 



In the case xq <S yo, ( |45| ) in the Appendix and (37) yield in the main logarithmic approximation 

X0\i I Uo \ , , J±xo 



E tot n-2w y [^\ log — + 4a,, — — . (39) 
\yoJ \a x J a x y 

We clearly see from the last expression that the energy is not bounded from below for xq ~ yo ~* 00 • The starting 
hypothesis xq <C yo is thus not satisfied. In the second case xq ^S> yo- With the help of ( |47| ) in the Appendix, the 
total energy per chain is 

E^t « -2wy log ( - 2w y ^\o g ( + 4a y ^. (40) 
\a x J x \ yo J a x y 

Minimizing with respect to xq and yo leads to xo ~ yo ~~ * oo, again in contradiction with the starting hypothesis 
x > y - 

These contradictions imply that the modulated line minimizes its energy at intermediate values of ro = Xo/yo, e.g. 
Tq ~ 1, where the Coulomb energy is of the order of the confinement energy. In this regime, (|37|) is difficult to treat 
analytically. 

We thus consider an ansatz giving the distribution of the kinks. Numerical simulations give zig-zag shapes. In this 
connexion we show in the Appendix that a sinusoidal profile 

y 

x(y) = xosin(27r — ) 

yo 

gives a correct description of the distribution in the continuum limit. The total energy is then given by 

4 ( l °l ai , 1 / 47rr |sin(t)| \ A J±_ 

Etot = -w y / dt K U J v 0= +4a y ^r , (41) 

n J°*/vo yt 2 + 167r 2 r2 S in 2 (t) \ J t 2 + Wir^ sm 2 (t) J a * 

where the upper cut-off corresponds to the dimensionless Debye length and K(k) is the complete Elliptic function of 
the first kind. The lower cut-off in the t— integral gives rise to the asymmetry between xo and yo, which is already 
present in (|^) and (flO]). As the sinusoidal distribution is valid in the continuum limit we shall take a x — > and 
dy — > keeping other lengths fixed. (^Tj) then depends only on ro . The numerical plot of this energy clearly shows a 
minimum for ro ~ 1.2, cf. Fig. 0, which agrees with the analytical predictions. 



9 




With the help of ( p6|) these results lead to 

xo ~ Vo ~ (w y /J±)a x > 1. 



(42) 



The condition that the lines do not overlap, xq <ti returns us to the boundary T <T C — \/v separating the phase 
of linear aggregates from the Wigner crystal of bikinks. 

A tilted stripe, which corresponds to a small argument of the sine distribution above, is easily shown to better minimize 
the energy. This fact is relevant in the context of cuprate oxides! 

WIGNER CRYSTAL OF SOLITONS 



ROUGH 
LINES 



WIONEK CRYSTAL OF B1KIN KS 



BIKINKS 

AND 

RODS 



LIQUID OF BIKINKS 



■0 



O C 
FIG. 3. Phase diagram of the 2D model. 
C is the confinement transition temperature below which solitons aggregate into bikinks. 0n is the crossover 
temperature below which bikinks aggregate into growing rods. At = the latter are roughened domain lines 
perpendicular to the chains when the Coulomb interaction r < T c . For < Q < On. this confined phase consists of 

aggregates screened by a liquid or a Wigner crystal of bikinks, depending on the strength of the Coulomb 
interaction. When T > r c both the lines at O = and the aggregates with their screening charges at < 9 < 0n 

crossover to a deconfincd phase, the Wigner crystal of solitons. 
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IV. NUMERICAL APPROACH 



In this section we consider the results obtained with the help of numerical methods. The calculations have been 
performed in the canonical ensemble with periodic boundary conditions. The latter implies interaction of each kink, 
within a computational cell, not only with all ether kinks but also with all their images residing in cells obtained by 
translation from the original computational celM. This results in an effective electrostatic interaction which is periodic 
in the computational cell used and which can be properly tabularized for rapid evaluationsEZI. At the beginning of each 
simulation we place the kinks at random and assign the sign to all of the Ising variables S n<a accordingly. We study the 
energy landscape of the system using the classical Monte Carlo technique with the standard Metropolis algorithm for 
the acceptance of kink motion. In order to rapidly reach configurations with the lowest energy, we perform simulated 
annealing from high temperatures with a variable sampling time range. We emphasize the common points between 
these numerics and the various points discussed above. Also, the impact of discreteness within the simulations is 
considered. 



A. The parameters of numerical simulations and the impact of the coarse graining 

The typical values taken in numerical simulations are J± w 10 -1 eV, v w 1% and a y ~ 4, 5A This leads to 
C w 100, Oo ~ 1, as above, and T sw 10/e. The Coulomb interaction is then monitored with the help of the dielectric 
constant. The value assigned to the coarse graining length needs some more care. It must be at least of the order 
of the core of a soliton (£o ~ 7 A in polyacetylene) . However, if it is too large the ground state of the system may 
be unreachable because of commensurability effects due to the discreteness of the lattice. The latter are especially 
relevant in the Wigner crystal phases where a pinning by the lattice occurs which distorts the Wigner crystal from 
its equilibrium position. We shall therefore estimate the upper boundary, allowing neglect of such commensurability 
effects, which can be reached by a x . 

We suppose that kinks form a Wigner crystal with their positions Xn 9 d in equilibrium, i.e. at the lowest energy 
distribution which is the triangular lattice&J. This structure melts when the fluctuations u of the kinks around their 
equilibrium position is a fraction 7, the Lindeman number, of their mean separation \J~sjv. The new positions are 

x n ,a — Xn% + u, with u = 7-y/ sjv , 7 < 1. Usually the melting is due to thermal or quantum fluctuations. However, 
if the coarse graining imposed by simulations is greater than u, kinks will also never reach their equilibrium position 
and the Wigner crystal will not be observed. Instead a two-dimensional commensurate phase will form. Thus, the 
coarse grained length must satisfy the condition a x < u — j-J £,QCi y /v. With the typical values given above this yields 
a x < 5O7A For this range of the coarse graining length the kinks will be able to adjust themselves and not be 
influenced by the underlying lattice. Supposing that 7 sw 0, 5 gives an overestimated value of the maximum length, 
a x « 25A 



B. Numerical results 

Figure |J shows the non- interacting phase at temperatures <C < ©o- This phase has been obtained from the 
interacting case by imposing a large dielectric constant e ~ 10 9 so that r -C 1. In accordance with what has been said 
concerning the non-interacting phase, a gas of bikinks is present and the latter aggregates into rods perpendicular to 
the chains. 
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FIG. 4. The non-interacting rod phase corresponding to <C B < Oo- 
This phase has been obtained by imposing a large dielectric constant. The computational cell displays the chains as 
horizontal lines, and their excitations as the points. Each excitation corresponds to the presence of a reversed spin 
with respect to the ferromagnetic ground state and thus represents a bikink. The aggregation of bikinks leads to the 

presence of rods perpendicular to the chains. 



Turning on the Coulomb interactions and adjusting the dielectric constant such that T 
a one dimensional Wigner crystal of domain lines as shown in figure |^. 



1 <C r c , we have at 6 = 
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FIG. 5. Domain lines at O = for a weak Coulomb interaction V ~ 1. 
As temperature approaches zero the rods of figure ^ cross the whole computational cell more like a system of 
domain lines. In the presence of a weak Coulomb interaction these lines order in the form of a one dimensional 

Wigner crystal. 



Turning o n stro nger Coulomb interactions the lines acquire a modulated shape in accordance with what has been said 
in section IIIC; the result is displayed in Fig. ^ for T w 10. By comparison with Fig. || the lines are still forming a 
periodic pattern. They roughen as a result of the competition between Coulomb interactions and confinement. 
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FIG. 6. Modulated domain lines at 9 = and K T < T c 
A Coulomb roughening of the domain lines takes place when Coul omb i nteraction is increased, in agreement with 

the results of paragraph [II C. 



V. CONCLUSION 



We have considered the effect of two competing long range interactions, Coulomb and confinement, on the thermo- 
dynamic properties of systems with two-fold degenerate ground-states. Analytical and numerical methods have been 
used and compared. The theory has been developed along the lines of the non-interacting case where a mapping to 
the Ising model allowed for an exact solution in 2D. Due to the complications brought by the Coulomb interactions, 
we have studied, in a semi-phenomenological way, the main features which distinguish the present case from the 
non-interacting one. In 2D the main parameter controlling the aggregation of solitons in the presence of the long 
range interactions is T in equation (pd]), i.e. the ratio between the Coulomb energy scale and the confinement scale. 
This parameter, below the crossover temperature 0o, cf. (^), is analogous to temperature above the confinement 
transition C , cf. (^J). The various confined and deconfined phases have been displayed on the phase diagram Fig. 
|3| depending on T and temperature. We have concentrated on the low temperature regime , i.e. below 0o, where 
we have shown that a crossover line T c ~ \jv separates a confined phase of bikinks and their aggregates from a 
deconfined phase which is a Wigner crystal of solitons. As expected, in the aggregated phases, e.g. rods and lines, 
the Coulomb interactions may play a crucial role because of the size dependent charging energy of the solitonic com- 
plexes. In the case of weak Coulomb interactions r -C 1, the efficiency of the bulk screening insures the preservation 
of the aggregates. The main effect of the long range interactions is to shift the chemical potential of the elementary 
defects by the Coulomb self-energy (p2[). However, when bulk screening is no longer efficient, we have found that 
antiferromagnetic interactions along the chains are necessary in order for the system to adjust the correct density of 
solitons and conserve the aggregates (p4|). The ground state of the system remains ferromagnetic to maintain the 
electroneutrality. In both cases, efficient and weak screening, a growing macroscopic aggregation takes place. Because 
of thermal roughening, only at T — do the complexes of solitons span the whole sample in the form of domain lines 
( |34| ) . We have shown that the latter are stable for larger Coulomb interactions 1 -C T < T c due to self-screening in the 
ensemble of solitons. This is manifested by the Coulomb roughening of the lines in accordance with numerical results, 
cf. Fig. ^|. The parameters of the distribution were derived from a variational procedure ([f2|). In the continuum limit 
a sinusoidal profile gives a correct description of the rough lines. 



APPENDIX 



We consider the periodic modulation for a general distribution of kinks x{y). The Coulomb potential energy associated 
with this distribution is given by (|38|). Going to the new variables t = (y — y')/2 and s — (y + y')/2 expression ( |3^ ) 
becomes 




1 



\J t 2 + {x{s + t) - x(s - t)f 



(43) 
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The integral has been regularized by subtracting the infinite contribution coming from the straight line x(y) = 0. 



In the case where xo <C yo, (43) can be divided as follows: 







rvo 




dt 1 


f ds 


J 4 


J + 


r 




1 a x 




Jxq 




a y | 



(44) 



In all integrals in ( fi~4| ) we can expand x(s ± t) « x(s) ± tx'(s). Thus the t and s integrals decouple and we obtain 



-2w„ 



log 



(45) 



In the opposite case, x$ 3> yo, we have 

2w 

f] , ~ 





/•xo 




dt J 


/ ds 


J ^ 






1 a x 









** l V /t 2 + (^(.S+t)-.T(.S-i)) 2 * 



(46) 



In the present case the expansion of x(s ± i) needs some care because xo 3> yo and the s and t integral do not easily 
decouple. The first t— integral, in (f46|), is over xo/yo 3> 1 periods. Its major contribution is at the vicinity of the zeros 
of x(s + 1) — x(s — t). These zeros exist for any s and correspond to "nesting" points t n = nyo/2 with n = 1, 2, ...2ro. 
From these considerations and with the help of 8t = t — t n , the major contribution to the integral reads 



2w y Y^ 



2r " ^dSt 



1 



% ^/tl + ASt^(s) \x'(s)\ ^ 



o 2r 



4x'(s) 



2w y / dnlog(n), 



where = xo/yo and we have taken x'(s) w ro. Thus, restricting ourselves to the main logarithmic approximation, 
we obtain 



-2w y log 



2^ — log 

xo 



In the second t— integral t is large and we can expand the denominator in x(s + i) — x(s — t). Going to the first order 
in l/ro the result is compensated by the regularizing term and thus gives zero. 
The Coulomb energy then reads 



E, 



coul 



yo 



Xo 



Xq 



yo 



2Wy lOg - 2Wy lOg 



•'•(I 



(47) 



The total energy of the system we are considering includes the potential energy and the confinement energy. The 
latter can be taken, for any xq, as 



j_l r\dx 

E C onf = a y — / — \dy = 4a y 
a x Jo dy 



Jj_xo_ 

a x yo 



(48) 



Combining (48) with either (45) or ([47|) the total energy per chain is thus given by (|37|). 



A particular function x(y) describing a single domain line in the two-dimensional plane at T = can be estimated. 
As implicitly assumed in (]48|), the confinement of this aggregate can be taken into account, in the Hamiltonian, with 
the help of a solid on solid model: 



H = — V - 6 +Vy"\x(y)-x(y-a v )\. (49) 

Minimizing (EjJ) with respect to x(y) yields the recurrence relation 
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~(y) = 



k 2 (sign (x(y) - x(y - a y )) - sign (x(y + a y ) - x(y))) + x (v') ( x (v) ~ x{y')f + (v ~ v') 



,/\ 2 



-3/2 



{x(y) - x{y')f + (y- y'Y 



-3/2 



(50) 



where sign is the sign function and k 2 — eV/e 2 = (tlodyT) -1 , has a unit of inverse length squared and is related to 
(|36|). The solution x(y) cannot easily be found analytically. However simulations give us a hint for they show a kind 
of zig-zag distribution. This suggests using a sinusoidal ansatz 



y 

x(y) = xosin(27r — ). 



(51) 



Comparing the resulting plot of each side of ( p0|) we can check the accuracy of this approximation. The result is 
displayed in Figure ([?]), which shows that the sinusoidal distribution is correct in the continuum limit. 



FIG. 7. The sinusoidal distribution. 
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